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We address the phase of a highly polarized Fermi gas across a narrow Feshbach resonance starting 
from the problem of a single down spin fermion immersed in a Fermi sea of up spins. Both polaron 
and pairing states are considered using the variational wave function approach, and we find that the 
polaron to pairing transition will take place at the BCS side of the resonance, strongly in contrast 
to a wide resonance where the transition is located at the BEC side. For pairing phase, we find 
out the critical strength of repulsive interaction between pairs above which the mixture of pairs 
and fermions will not phase separate. Therefore, nearby a narrow resonance, it is quite likely that 
magnetism can coexist with s-wave BCS superfluidity at large Zeeman field, which is a remarkable 
property absent in conventional BCS superconductors (or fermion pair superfluids). 



Whether an s-wave superconductor (or fermion pair 
superfluid) can coexist with magnetism is a long-standing 
issue in condensed matter physics. Back to 1960s', Chan- 
drasekhar and Clogston independently considered the re- 
sponse of a BCS superconductor to spin polarization due 
to a Zeeman field [U [2] ■ They found that an s-wave su- 
perconductor will remain unpolarized until a critical Zee- 
man energy h c ~ A/y2 with A being the pairing gap, at 
which the system undergoes a sharp phase transition to 
a partially polarized normal state, and this critical field 
is now known as Chandrasekhar-Clogston (CC) limit of 
superconductor. In this scenario, superconductivity can 
not coexist with magnetism. Later on, there are several 
proposals for magnetized s-wave superconducting states, 
and the most famous ones are the Fulde-Ferrell-Larkin- 
Ovchinnikov state [3] and the Sarma State [3]. However, 
so far none of them have been firmly observed in an s- 
wave BCS superconductors. 

In the last few years, this problem has been revis- 
ited by a series of experiments on two-component Fermi 
gases with population imbalance [5-9 . Experiments have 
reached a consensus that at resonance regime and the 
BCS side, there is a CC limit, where a direct transition 
from a fully paired fermion superfluid to a partially polar- 
ized normal state has been observed, and no evidence of a 
magnetized superfluid has been found [51 18141 0|. However, 
all these studies were done across a wide resonance. Re- 
cently, several experimental groups have begun to study 
narrow Feshbach resonances, such as 6 Li at 543. 25G [TT] 
or 6 Li- 40 K mixture at 154.719G [T2], where the effect of 
a finite resonance width needs to be taken into account. 
In this letter we find that the resonance width indeed 
has dramatic effects on physics of highly polarized Fermi 
gases. 

In contrast to a wide resonance, to character a Fermi 
gas nearby a narrow resonance one not only needs 
hpa®! where a® is zero-energy scattering length between 
fermions, but one also needs to consider h 2 ftp / (2mW a^g) 
with m being atom mass W the resonance width and db g 
the background scattering length. If h 2 kF/(2mWah g ) > 
1 the resonance is considered to be a narrow one. Another 
important quantity is kpa^b with abb being the scatter- 
ing length between closed channel molecules, the effect 



of this quantity will be analyzed later. In this work we 
focus on highly polarized limit and show all above three 
parameters play an important role in determining the 
nature of many-body phases. The studies of this work 
contain two parts: 

Firstly, we consider a single down spin immersed in a 
Fermi sea of up spins. Two different types of states are 
compared, which are polaronic state and pairing state. 
For polaronic state, the single down spin is dressed by 
particle-hole pairs of up spins, and becomes a fermionic 
quasi-particle [F3] . If this state has lower energy, the sys- 
tem will be a normal state of polaron Fermi liquid at 
sufficient high polarization. For pairing state, one of the 
up spins will form a bound state with the single down 
spin. If this state has lower energy, each down spin will 
form a pair, and the system will be a mixture of con- 
densed pairs and majority fermions. For wide resonance, 
a polaron to pairing transition takes place at the BEC 
side of the resonance [T4"Ml6j . Here we show that as the 
width of resonance gets narrower, the transition point 
will be shifted toward the BCS side. We find out how 
the critical value of (fcpeig) c changes with the quantity 
h 2 k F /(2mWa hg ). 

Secondly, when the pairing state has lower energy, the 
mixture of pairs and fermions may phase separate due 
to the repulsion between pairs and fermions. A strong 
enough repulsion between pairs are crucial to stabilize a 
uniform mixture. For a given H 2 fcp / (2mW a^ g ) , we find 
out the critical repulsion (fcpdbb)c as a function of h-pa®- 

Hence, we conclude that when l/(kpa^) > l/(kpa^) c 
and fcpabb > (fcFflbb)c, it is energetically favorable for mi- 
nority fermions to form pairs, and condensate of fermion 
pairs can uniformly mix with majority fermions, that is 
to say, magnetism can coexist with fermion pair super- 
fluids at highly polarized Fermi gases. The fact that this 
can happen at the BCS side and resonance regime repre- 
sents a significant distinction between narrow and wide 
resonances. As far as the response to spin polarization 
is concerned, at resonance, or even at the BCS side, this 
system behaves similar as the BEC side of a wide reso- 
nance. This picture is also consistent with a recent high- 
temperature study of narrow resonance [T7] . 

Model: We use the following two-channel model H = 
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Hq + V c + Vb g + Vbb to describe a narrow resonance 
Ho = $>k + "o)btb k + + 44^) (1) 

k k 

V c = g X] A k(6qrfq-k«k + Wkdq_ k 6q) (2) 
kq 

Vbg = U ^ A k A k /ttj ( .(i q _ k d q _ k /Uk/ (3) 

kk'q 

Vbb = ^^2 ^ q _ k &q-k'6k' (4) 
kk'q 

where v) and d) are creation operators for majority up 
spin and minority down spin, respectively, and is 
creation operator for bosonic closed channel molecule. 
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7 = m /m u is the mass ratio. V c and Vb g represent 
the interchannel coupling and the background scatter- 
ing, respectively. A k = 0(A — |k|) and A is the mo- 
mentum cutoff. In the Hamiltonian, the molecule de- 
tuning i>q, the inter-channel coupling g and the back- 
ground interaction parameter Uq are bare quantities with 
A dependence, which need to be renormalized as follows 
"b(A) =vr-[l- Z{A)}g 2 /U^ 50 (A) = Z(A)g T , 
and U Q {A) = Z(A)U Z , where Z(A) = (1 - U^A/tt 2 )- 1 , 
l/m r = l/m u + l/m d and U r — 2 , Ka^ & /m T . The renor- 
malized quantities U r , g T and v r are related to a s as 
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and the zero-energy scattering length a° is given by a° = 
m T (U r — g 2 1 v T ) I (2tt) . Denoting v r = Afi(B — B ), where 
B is the location of the resonance and A/j, is the differ- 
ence of magnetic moment between two channels, and in- 
troducing W — g 2 /U r , we have a° = a h<r {l — W/[Afi(B — 
B )}}, and a s (E) = a bg {l + W/[E - A(i(B - B )}} 
where Wab g is always positive. For C/bb> since we only 
consider it to the mean-field order, we will take it as 
= Airh 2 db b / '(rn u + m d ). Our following results will be 
presented in terms of physical parameters (W, Bq, fcpBg, 
^FObg and fcpObb)- 

Polaronic State: We first adopt the following varia- 
tional wave function which includes one particle-hole con- 
tribution 
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Here and below, all the summations with the prime (') 
of k and q are restricted to |k| > fcp and |q| < fcp, 
respectively. After energy minimization we obtain a self- 
consistent equation for polaron energy 
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FIG. 1: (Color online) Polaron energy as a function of 
fi 2 fep/(2m"W / ab g ). (a) 7 = 1 and different interaction pa- 
rameters. l/(kpas) = (dash line) and l/(kFa°) = —2 (solid 
line), (b) l/(fcFOs) = with different mass ratio. 7 = 6/40 
(solid line), 7 = 40/6 (dash line) and 7 = 00 (dash-dotted 
line). All curves are results from one particle- hole approx- 
imation, while solid black dots show results including two 
particle-hole pairs contributions. fcFCtbg is chosen as —0.1. We 
will also set h = 1 in all figures to simplify the presentation. 



where r 2 (q, E + e q ) coincides with two-particle vertex 
with total momentum q and total energy E + e q within 
the ladder approximation. This is because the variational 
wave function \ip p } describes the processes that an up 
spin (taken out from an occupied state q) and the single 
down spin undergo repeated scattering, as well as coher- 
ent conversion between open and closed channels. This 
physical process is precisely what is captured by the lad- 
der approximation. 

The explicit form of r 2 (q, E + e q ) is given as 



r 2 ( q ,i;+e q ) = 
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where I(q,E + e q ) = £ k l/[c£ + e q _ k - (E + e q )] - 
E k V( e k + e - k )- The difference between Eq. <TrT) and 



the previous results for a wide resonance [T31 Q3] is that 
a constant a s is replaced by an energy dependent one 



a s {E- 



where E + e q — e q represents the energy 



of the relative motion for two atoms undergoing repeated 
scatterings. 

The polaron energy as a function of h 2 kF / (2m u Wab g ) 
is plotted in Fig. from solving self-consistency equa- 
tion Eq. m. As one can see, when the dimensionless 
parameter fi^k-p / {2m u W a^g) increases from zero, (i) for 
7 = 1, and nearby resonance l/{kpai) ps polaron en- 
ergy E will increase, while at the BCS side l/(kpa^) <C 0, 
E will decrease; and (ii) at resonance, if 7 is greater 
than a critical value, E will also decrease. We have 
also checked the energy convergence by considering two 
particle-hole contributions [2QJ [2Jj| ■ The numerical solu- 
tions with two particle-hole contributions are shown as 
the dots in Fig. ([!]). For 7 = 1, one can see that the 
corrections from two particle-hole pairs are always negli- 
gibly small, and it becomes even smaller as \ W\ decreases. 
While for 7 — > 00, the deviation is a little larger (dashed 
line and the dots in (b)), as already noted for a wide res- 
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onance in 20J, but it is still within only a few percents. 
This result justifies the validity of the expansion in terms 
of the number of particle-hole pairs in computing energy 
for a narrow resonance. 

Pairing State: For pairing state, we use the variational 
wave function first introduced in Ref. 1141: 
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where |FS ) refers to the Fermi sea with one spin-f parti- 
cle removed from the Fermi surface of |FS) in the polaron 
state (JsJ . If we only consider bare pair wave function 
without including particle-hole contribution, the pairing 
state energy is given by r^" 1 (0, E + 6f) = 0, as shown in 
the lines of Fig. The interaction between pair and 

majority up spins can be described by including particle- 
hole pairs. Up to one particle-hole pair, by minimizing 
energy we obtain a closed integral equation [3T] , and the 
numerical solution of these equations are also shown in 
Fig. ([2]). We find that as kp / '(2m u W ab g ) increases, the 
pairing state energy decreases at the BCS side and at 
resonance regime, despite of different mass ratios; while 
it increases at the BEC side. Another important feature 
one can find from Fig. ^ is that in the limit W — > 0, the 
pairing state energy always saturates to —ep. This can 
be understood as follows: when one down spin is added 
into the system, an up spin is taken out from the Fermi 
sea (subtract energy ep) to form a pair with the down 
spin, whose energy approaches v t in the limit W — > 0. 
Thus the pairing state energy should approach — ep + is r 
|24j . Moreover, at any fixed a° the ratio W/v T is fixed, 
thus v T — ¥ and therefore pairing state energy always ap- 
proaches — tp as W — > 0. This also indicates that the in- 
teraction between a pair and the residual majority atoms 
vanishes in the limit W — > 0. We have performed a three- 
body calculation and find out the atom-dimer scattering 
length a a b from the asymptotic behavior of the three- 
body wave function [5TJ [25] , as shown in the inset of Fig. 
([2]) (a), which indeed shows a a b — > as W — > 0. 

Polaron- Pairng Transition: The transition from pola- 
ronic state to pairing state can now be determined by 
comparing their energies. In Fig. (|3])(a-b) we consider 
two concrete samples studied in current experiments: 6 Li 
at 543.25G and 6 Li- 40 K mixture at 154.719 G, and the 
parameters are typical values taken from experimental 
papers [TTJ Q21 [21] ■ We found that in both cases, the 
polaron to pairing transition is located at the BCS side 
of the resonance, which is away from resonance with 
A/j,(B — Bq) on the order of ep. At the transition points, 
l/(fc F a°) = -4.35 for 6 Li and l/(k F a°) = -0.55 for 6 Li- 
40 K mixture, where the systems are very BCS-likc. This 
transition has also been observed in a recent experiment 
on 6 Li- 40 K mixture and the transition is indeed observed 
at the BCS side [22]. 
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FIG. 2: (Color online) Pairing state energy as a function of 
H 2 k F /{2m u Wa bg ). (a) mass ratio 7 = 1 and different inter- 
action parameters. l/(fcpas) — (dash line), l/(fcpa!?) = 
— 1 (solid line) and l/(kFa®) = 1 (dash-dotted line), (b) 
1/(/df<1s) — but different mass ratio. 7 = 6/40 (solid line), 
7 = 40/6 (dash line) and 7 = 00 (dash-dotted line). All the 
curves are computed from bare pairs, while the dots, crosses 
and triangles are results with one particle-hole contributions. 
Inset of (a): the atom-dimer scattering length a a b (in unit of 
\/m u Eh/h with E\, the two-body binding energy) as a func- 
tion of v r /W. fcpabg is chosen as —0.1. 



In Fig. |3][c-d), we plot the critical value 
l/(fcpfls) f° r polaron to pairing transition as a func- 
tion of kF/(2m u Wab g )- One finds that when 
h 2 k-F/(2m u Wab g ) > 1, the transition will be shifted to 
the BCS side. This condition is equivalent to ^5 
1/lfcpabgl- Since usually fep^bg <C 1, it means that the 
resonance width does not need to be very narrow. One 
also notes that the transition point is not sensitive to 
the value of /cFabKjtself (Fig. (|3])(c)), but is sensitive to 
mass r atio (Fig. jSjd)). The inset of Fig. ^(c) shows 
that in the limit W — > 0, the critical point will approach 

— > ep, which means that the pairing state will be fa- 
vored once the energy of closed channel molecule is below 
the Fermi energy. 

Stability of the Mixture: Discussion above concludes 
that for sufficiently narrow resonance, the highly polar- 
ized Fermi gases contain a mixture of majority fermions 
and bosonic pairs. Next question is whether they will uni- 
formly mix or phase separate. To answer this question, 
we note that for very low density of down spins and suf- 
ficiently narrow resonance, we can expand the cquation- 
of-state in terms of the density of bosonic pairs rib up to 
the second order 261 



(9) 



where £p is the energy density of spin-f Fermi sea, 
/ijj = E+e p where E is the pairing state energy computed 



above. The repulsion g = m t^ m d Qbb + ffmd contains 
the contribution from the bare interaction between closed 
channel molecules and the induced interaction g- m ^ from 
the inter-channel coupling, which is calculated within 



4 



-0.2 



LJJ 



-0.4 



-0.6 




0.35 0.7 
-1/(k F a ^ 



-0.5 



-1 



(b) 



2.5 




k^/(2m u Wa l , 
F v bg 



FIG. 3: (Color online) Comparison of polaron energy (solid 
line) and pairing state energy (dash line) as a function of zero- 
energy scattering length for 6 Li- 40 K mixture (a) and 6 Li (b). 
For (a), we take k F a bg = 0.022 and W = 54.91e F [22], and 
for (b) we take k F a hg = 0.016 and W = 12.2e F [11] ■ In 
(a), 40 K is taken as minority component. Panels (c) and (d) 
show critical value of the transition l/(k F a°) c as a function 
of h 2 k F /(2m u Wab g ). For (c), 7 = 1 but different fcpabg. 
fcpabg = —0.1 (solid line) and k F a bg = 0.1 (dashed line). For 
(d), fcFGtbg = 0.01 but different 7. 7 = 6/40 (solid line), 7 = 1 
(dash-dotted line) and 7 = 40/6 (dash line). Inset of (c), the 
critical value of transition in term of (in unit of tp) as a 
function of h 2 k F /(2m u Wa bg ). 



The stability condition against phase separation is given 
" > EH, from whieh we ean de- 



by 



tcrmine the critical value for abb- The results are plot- 
ted in Fig. [i] for 6 Li- 40 K mixture (a) and 6 Li (b) in 
the regime where pairing state is favorable. We can see 
that it requires fcpObb > 0-81 for 6 Li- 40 K mixture, and 
fcpabb > 0.017 for 6 Li at resonance. Very likely, this con- 
dition can be satisfied in 6 Li but not in 6 Li- 40 K mixture. 
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Note Added: Recently, we became aware of the ex- 
perimental work from Innsbruck group in which the po- 
laron properties and the polaron to pairing transition has 

been observed [221 ■ We also became aware of two other 
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Born approximation [21]. From Eq. ^ we obtain: 
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FIG. 4: (Color online) Critical value required for k F a bb to 
prevent phase separation as a function of l/{k F di), for 6 Li- 
40 K mixture (a) and 6 Li (b). 
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I. SUPPLEMENTARY MATERIAL 

A. Polaron energy including two particle-hole contributions 

In this section of supplementary material, we derive the self-consistent equation for polaron energy including two 
particle-hole excitations. 

The Hamiltonian under consideration is given as equation (l)-(3) in our paper, and the variational wave function 
is written as: 
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where \FS) is the Fermi sea of N spin-f particles, |0 k q) and |?y q ) correspond to single particle-hole excitation while 
| a), |/3) represents the contribution from two particle-hole excitation. All momentum summations are restricted 
within |q| < fcp, |k| > as mentioned in our paper. Now we solve the Schrodinger equation H\ip) — E\ip) and drop 
all the terms beyond two particle-hole excitation, we get the following set of equations: 
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where </? q = U J2'k ^kq, G kqq - = AU J2' K a k K qq ' and we have defined: 
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Finally, after eliminating all the variables other than Gk qq ' from equation ( 15 l-([T9|) and using the renormalization 
relations between Uq, go, vq and U r , g r , v r as given in our paper, we obtain the following integral equations: 



kq' 



Tkqq'Gkqq' — </>kq' — </»kq + E Gk'qq' /-Ekk'. 



»kq : 



E, 



kq 



q' V k'q' q 



(20) 
(21) 

(22) 



where 7 k qq ' - ^(q + q' - k, E + e q + e», - e k ), £ kq = e k + <-k " % ~ E and £ k k' qq ' = ^ + e«, + e q+q ,_ k -k< - 
£q — e q / — -E. r^ 1 has been defined in equation (6) in our paper. The above integral equations are solved by numeric 
iterations to obtain the polaron energy E. 



B. Three-body problem and the atom-dimer scattering length 

In this section of supplementary material we obtain the atom-molecule scattering length by solving the three-body 
problem using a two-channel model (see equation (l)-(3) in our paper). We first write he three-body wave function 
in the following second quantized form in momentum space: 



1-03 



E < Wa u k 1 u k 2 d -k 1 -k 2 +J2 T 'P U l b - 
Lkik 2 P 



| vac) 



From the Schrodinger equation H\ip 3 ) — E\-ip 3 ), we have: 

{el + el, + e£ +k , - E)<f>w + U E^pk' " ^ P k) + \go(w - %) = 

p 

{el + 4 + v - E) Vk + 2g E 0k'k = 



(23) 

(24) 
(25) 



Defining — Uq </>k'k and using the renormalization relation (here the introducing of the prefactor Uq is just 
to eliminate the ultra- violate divergence in 53k' > this divergence corresponds to the - singularity in the real space 
wave function), we have: 



7kC*k 



E 



Ok' 



k'k 



>/k 



where E k = E - el - e\ - v , E k , k = el + eg, + e£ +k , 

9 2 



7k = [U, 



2ffr «k 
U r Ek 

E and we define 

-l 



E-el-ei-v r 



E 



i i 



^kk' e r k , 



(26) 
(27) 

(28) 



The fourier transform a(r) = {2ir) 3 J d 3 kake lk r is just the asymptotic wave function between an atom and a molecule. 
To see this, we first write wave function (23 1 in real space: 



|V>3&) = J d 3 ri d 3 r 2 d 3 r 3 [0 o (r 1 ,r 2 ,r 3 )u t (r 1 )ii t (r 2 ) ( i t (r3) + ( /) c (r 1 ,r 2 ,r3)u t (r 1 )6 t (r3)] |vac) (29) 

Mri,r 2 ,r 3 ) = E ^k^^)^^) (30) 

kik 2 

* c (ri,r 2 ,r 3 ) = 6 3 {r 2 - r 3 ) E^e^"" 1 "^ = 6 3 {r 2 - - r 3 ) (31) 
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FIG. 5: Induced molecule-molecule interaction under Born approximation. The dashed and solid line refers to the propagator 
of molecule and atoms respectively. Each vertex is given by the inter-channel coupling g r . 



In the limit r 2 - r 3 -> 0, we have o (ri,r 2 ,r 3 ) ~ £ ki (S k , 0k 1 k 2 )e' kl (ri r3) = jj^a(ri - r 3 ) and (j) c (ri, r 2 , r 3 ) = 
5 3 (r 2 — r 3 )?7(ri — r 3 ). On the other hand, from equation ( |27| ) we can see that 7?k and ctk has identical singularity as 
k — » (since -J- is analytical near fc = 0). This means that in real space, ry(r) and a(r) have the same long range 
behavior as r — > oo. As a result, the long range behavior of ev(r) represents the asymptotic wave function between an 
atom and a molecule as they separate far away from each other. 

As r — > oo, we have a(r) = 1— s ^r IL (l + 0(7' _1 )) where a am is the atom-molecule scattering length. This leads to the 

fact that ak = <5 3 (k) — 5^2 2 p E -(l + 0(fc)) as k — > 0. As a result, it will be very convenient to define «k = S 3 (k) 
and we have: 



1 fik 
2tt 2 A; 2 



/3k = 



-fc 2 

47T7k 



#ok 7 ^q 2 E^ 



(32) 



The above integral equation is solved numerically and the atom-molecule scattering length is given as a a 



C. Induced Interaction 



In this section of supplementary material, we evaluate the induced interaction g m< i within Born approximation. 
The corresponding diagram for the T-matrix of molecule-molecule scattering is shown in Fig. [5] We have <?; n d = 
4 ^l ai " d = T m and: 

_ 4 f dw f ri 3 p 1 1 

^-f-y 27r y (27r)3 e t )2(jw+e i )2 w 

where g 2 = 27r ^ gM/ , = ^ an d M4. = is given by the molecular energy in the single impurity problem. 

All the external frequency and momentums are set to zero. 

After performing the frequency integral and some change of variables, we finally have: 

k F a ind = - 7 (1 + y-^fkrobg—) V( 7 , i^l) (34) 
7r V " ep / ep 

f°° x 2 1 

F ^ = 2 J 1 ^[(l + 7 -i), 2+ ,-ip-2( 7 - 1+ ,) 2 (35) 

where 7 = — . 



